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ABSTRACT 

Knowing the site of 7-ray emission in AGN jets will do much for our understanding of the 
physics of the source. In particular, if the emission region is close to the black hole then absorption 
of 7-rays with photons from the broad-line region could become significant. Such absorption is 
predicted to produce two specific spectral breaks in the 7-ray spectra of Flat Spectrum Radio 
Quasars (FSRQs). We test this hypothesis using 3 years of Fermi observations of nine bright 
FSRQs. A simple power law fit to the spectrum of each source can be significantly improved by 
introducing a break, but the break energies are inconsistent with those predicted by the double- 
absorber model. In some cases the fit can be further improved by a log-parabola. In addition, 
by dividing the data from each source into two equal epochs we find that the best description of 
an object's spectrum often varies between a log-parabola and a broken power law. 

Subject headings: BL Lacertae objects: general; galaxies: active; galaxies: jets; gamma rays: general; 
radiation mechanisms: non-thermal 



1. INTRODUCTION 

The commonly accepted paradigm is that ac- 
tive galactic nuclei (AGN) consist of a central 
black hole surrounded by a torus of dust. There 
are also two regions of material called the broad- 
line region (BLR), which extends to a distance of 
< 0.1 pc from the cent ral black hole for bright 
AGN (iKaspj et al.lll996l) . and the narrow-line re- 
gion (NLR), which is loca ted ~ 10 - 100 pc from 
the central black hole ( Alexander fc Hickox,2012,) . 
Some AGN also have a relativistic jet perpendicu- 
lar to the disk. If this jet has a small line of sight 
to the observer the object is termed a blazar. 

Blazars emit radiation across the whole electro- 
magnetic spectrum including the 7-ray regime; for 
example, the Fermi Gamma-ray Space Telescope 
has detected hundreds of blazars at photon ener 



Fi, is the flux at frequency v. Blazars can be 
divided into two subclasses. Flat Spectrum Ra- 
dio Quasars (FSRQs) and BL Lacertae objects 
(BL Lacs). BL Lacs tend to have their sec- 
ond peak at higher energies than FSRQs which 
tend to have a secon d peak below ~ 100 MeV 
( Ghisellini et al. I I2OO9I) . [n BL Lacs emission lines 
from the BLR tend to be weak o r not-observed 
( Georganopoulos fc Marscher 19981 ). 

The 7-ray emission from blazars is assumed 
to originate in the jet, and observed variabil- 
ity timescales suggest that the 7-r ay emission re- 



gion is small in size (< 10 ^ pc) (jSbarrato et al 



20111 ). However, both the region and the mech- 



gies be tween 100 MeV and 100 GeV (jAbdo et al 
2010b( ). The shape of a typical blazar's spectral 



anism of 7-ray emission is currently unidenti- 
fied. One popular theory is that electrons in the 
jet emit seed photons via synchrotron emission 
and inverse Compton scatter these photons up 
to 7-ray energie s (the synchrotron se lf-Gompton 



energy distribution is characterized by two peaks 
when plotted in units of log(j^) against \og{uYy) 
where v is the frequency of the radiation and 



model, see e.g. iMaraschi et al.l Il992f) . There is 



also widely thought to be a contribution from pho- 
tons that originate outside the jet that become up- 
scattered (the external-Gompton model, see e.g. 
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Dermer fc Schlickeiseij|l993l ). This contribution is 
thoug ht to be more importa nt in FSRQs than BL 
Lacs (jGhisenini et alJll998l ) due to the availabil- 
ity of photons from the BLR for upscattering. In 
these models the synchrotron emission accounts 
for the first, low energy peak in the SED while the 
Compton emission accounts for the second, high 
energy peak in the SED. Alternatively, hadronic 
components in the jet may be accelerated to high 
enough energies to emit synchrotron radiation in 
the 7-ray regime which accounts for the high en- 
ergy peak in the SEP (hadronic models, see e.g. 
Miicke et al.ll2003l) This is still an active area of 
research (see e.g. Bottcheij I2OIOI for a short re- 
view) since it is key to extending our knowledge of 
the underlying jet physics. In order to construct 
and test these models it is crucial to understand 
the 7-ray emission region and the shape of the 7- 
ray energy spectrum being produced. 



Abdo et al.l (|2010al ) concluded that the spectra 



of FSRQs differ from simple power law descrip- 
tions in the Fermi energy range (~ 100 McV - 
300 GeV). In particular, they concluded that for 
the brightest FSRQ source, 3C 454.3, a broken 
power law gives the most acceptable fit compared 
with si mple power law an d curve d (log-parabola) 
forms. Poutanen fc Stern ( 2010l hereafter PSIO) 
note that recombination of electrons with nuclei 
can produce sharp increases in photon density 
above a specific energy. Therefore the opacity to 
7-rays due to photon-photon (7-7) pair produc- 
tion would also see a sharp increase above a cor- 
responding energy. PSIO used a simple photoion- 
ization model of the BLR clouds and an emission 
spectrum i n the energy range 1 eV - 100 keV from 
Laor et al. 1 ([1997) to model the opacity the BLR 
presents to a 7-ray of given energy. They con- 
cluded that there would be two increases in opac- 
ity, one occurring in the region 4 GeV - 7 GeV 
due to He II recombination and one in the re- 
gion 19.2 GeV - 30 GeV due to H I recombina- 
tion. A further simplification in PSIO predicted 
that these changes should appear at 4.8 GeV and 
19.2 GeV respectively. It was noted that if a sim- 
ple power law of 7-rays is subjected to such an in- 
crease in opacity then the effect is a change in the 
spectral index akin to a broken power law. Fur- 
thermore, it is suggested that this is the cause of 



the bro ken power law spectra found in I Abdo et al 



higher energy. This scenario is termed the 'double- 
absorber model'. In PSIO, the authors favor the 
double-absorber model over fitting a power law 
with a single break in which the break energy is 
left as a free parameter. This is a very impor- 
tant result if correct, since absorption of 7-rays by 
BLR photons strongly implies a 7-ray production 
region that lies within the BLR. 

We perform analysis of the Fermi data from 
several bright FSRQs in order to test three spec- 
tral forms to determine which gives the best de- 
scription of the data. The first is a simple power 
law 



dE \Eo 



(1) 



where dN /dE is the differential photon flux as a 
function of photon energy, fc is a normalization 
constant, E^ is a normalization energy, which we 
fixed at 0.5 GeV, and P is the spectral index. The 
second form tested is a broken power law, which 
describes a sudden change in the spectral index of 
a power law from Pi to P2 at a break energy Ei, 



dE 



Er, 



(2) 



if E>Eb 



If notable deviation from a simple power law was 
found to be present, we tested to see if the break 
energies were in agreement with those predicted 
from the double-absorber model and if they were 
stable over time. The final form tested is a log- 
parabola, a smoothly varying function that natu- 
rally also describes a peak frequency 

-r-/3iog(J-) 



dN 
'dE 



(3) 



(|2010af ) and that a second break should occur at a 



where we again fix the normalization energy to 
0.5 GeV and /3 is a parameter that defines the 
curvature of the spectrum. A power law is phys- 
ically untenable over a sufficiently wide range of 
energies whereas the log-parabola can naturally 
describe the curvature expected on either side of 
the peak in the SED. 

We used a larger dataset than was available 
for the PSIO analysis and also made several im- 
provements to the methods used in PSIO. As the 
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authors in that paper noted, PSIO used a sim- 
phfication in the Fermi-LAT's containment angle 
and effective area which degraded the accuracy 
of photon selection. Secondly, PSIO used only a 
single region of space to characterize the energy 
spectrum of the 7-ray background for each ob- 
ject rather than a more thorough model of diffuse 
7-ray emission and emission from known nearby 
point sources. Thirdly, PSIO used a binned anal- 
ysis which introduced extra degrees of freedom in 
the placement of the bin edges. These simplifica- 
tions had a noticeable impact on the precision of 
the fits. In the broken power law fits performed in 
PSIO the uncertainties in the break energy were 
on average 50% larger than those for the same 
objects ove r the sa me time period performed in 



Abdo et all (|2010al ). 



Our improvements to this analysis are detailed 
in Section 2 along with the selection cuts we 
applied to the data, how we tested for breaks 
and curvature in the spectrum, and whether any 
breaks were stable over time. Our results are pre- 
sented in Section 3 and our conclusions are pre- 
sented in Section 4. 

2. Data Analysis 
2.1. Event Selection 

We analyzed the 9 FSRQs which were selected 
in PSIO for their brightness and lack of contam- 
ination from other point sources or the Galactic 
bulge. These sources are listed in Table 1. 

We examined all photons that arrived within a 
15° radius of the target using the P6_V3.DIFFUSE 
version of FermVs response function. Following 
standard practical], we required that a photon's 
energy and incident angle on the LAT were suit- 
ably well reconstructed (event class 3 (diffuse) or 
better) and any photons received when Fermi's 
zenith angle was < 105° were discarded to avoid 
7— ray contamination from the Earth. 



ergy range where a spectrum is well described by 
a simple power law except for a single break (be 
it caused by absorption or otherwise) then it is 
straightforward to identify this break. However, 
consider the case where the intrinsic spectrum (i.e. 
neglecting the break) is so curved that it is not well 
described by a simple power law. In this case, the 
break energy identified by fitting a broken power 
law would try to describe both the intrinsic cur- 
vature and the true break. In this way the break 
energy found by the fit is biased from the true 
break. Because of this it is important to identify 
an energy range where curvature in the spectrum 
is not too great. We performed Monte Carlo sim- 
ulations using the standard tool gtohssim of an in- 
trinsically curved spectrum with a spectral peak at 
^ 250 MeV (the spectral parameters were chosen 
to match the greatest curvature found in the obser- 
vational sample, F = 2.11 and /? — 0.12) and intro- 
duced an absorption feature above 1.8 GeV that a 
broken power law fit should correctly identify (an 
optical depth tt = 2 was used which corresponds 
to AF = 0.5 in a spectrum with no curvature). 
When we chose the spectral peak as the threshold 
of the energy range, the break energy was cor- 
rectly identified (1.7 ± 0.2 GeV). However, if we 
reduced the energy threshold below the spectral 
peak, the break energy was misidentified (in half 
the simulations, the break energy varied by more 
than 1 standard error if the energy threshold was 
reduced in 50 MeV steps). When we removed the 
absorption feature from the spectrum, a thresh- 
old of 50 - 100 MeV above the peak was some- 
times needed to find a stable break. The break 
energy found varied widely between simulations 
(3 ± 1 GeV) demonstrating degeneracy in fitting 
a broken power law to a curved spectrum. These 
simulations show that if there are spectral breaks 
in a curved spectrum they can be correctly iden- 
tified provided the threshold energy is above the 
spectral peak. The spectral peaks of each object 
i n our sample have been estimated by lAbdo etaT 



2.2. Testing For Deviation from a Simple 
Power Law 

Close to the high energy peak in the SED, we 
expect there to be significant curvature in the 
spectrum of any given object. If we have an en- 



(j2010cf) from the multiwavelength SED. However, 
there is some uncertainty as with any observation 
and it is also possible that the spectral peak has 
shifted with time, meaning the actual peak could 
be higher than these estimated values. With this 
in mind, for each object in our sample, we fit- 
ted a broken power law with the threshold first 

^ |http://fermi.gsfc. nasa.gov/ssc/data/analysis/scitools/likelihoodItutaiti atikanj ieak energy and then 50 MeV higher to 
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check the break was stable. If it was not, we in- 
creased the threshold in 50 MeV steps until a sta- 
ble break energy was found. The minimum energy 
determined in this way was used when fitting all 
spectral forms. In all cases the threshold was well 
below where the double-absorber model predicted 
breaks to occur so they would be identifiable by 
fitting a broken power law over the appropriate en- 
ergy ranges (see the following section). Since only 
bright sources are investigated in this work, the 
loss of photons incurred by increasing the energy 
threshold is not detrimental, especially since the 
containment angle of photons is relatively poor at 
the low end of the spectrum meaning the signal-to- 
noise ratio is degraded. In the case of PKS 2022-07 
we did not use the peak energy from Abdo et all 
()2010c^ since it could not be estimated directly 
from the object's SED but could only be estimated 
from the empirical relationship between spectral 
index and SED 



logKcak) = -4.0-r + 3.16. 



(4) 



For this object we fitted all of the data above 
100 MeV with a log-parabola in order to determine 
the peak energy, which we found to be 290 MeV. 

For each object we took all the photons pass- 
ing the selection cuts from the low energy thresh- 
old up to 100 GeV. We then used gtlike to de- 
termine the likelihood of the data being drawn 
from the three spectral models: simple power law, 
broken power law and log-parabola. Because of 
the low energy threshold imposed, break energies 
close to this value could not be tested for. The 
likelihood calculation incorporated point sources 
in the region of interest with parame ters fixed to 



those in the Fermi 1 Year Catalog (jAbdo et al 



2010br ). Galactic and extragalactic emission were 
accounted for using the models glLiem_v02 and 
isotropic_icm_v02 respectivel}!!. The Akaike In- 
formation Criterion (AIC) of each spectral form 
was then calculated from its likelihood. The dif- 
ference in AIC between models estimates their rel- 
ative Kullback-Leibler information quantities, i.e. 
how much each model tested diverges from the 
distribution of the data relative to the other mod- 



ances on the one hand the random error in the esti- 
mator from overfitting which results from a model 
containing more free parameters, and on the other 
hand the systematic error in the estimator from 
the bias of a mo del containing fewer parameters 
(|Bozdoganlll"987h . The spectral form with the low- 
est AIC value is taken to be the best description of 
the data and differences of 2 or more between AIC 



value s are considered significant (see iLewis et al 



201 iL and references therein). If a simple power 
law was not an adequate fit we tested if the break 
energies were in agreement with those predicted 
by the double-absorber model and if the break en- 
ergies were stable with time. 

2.3. Testing For Energy Breaks Caused By 
BLR Pair-Production 

As noted earlier, it is pointed out in PSIO that 
a photon-photon absorption feature in a power 
law spectrum can be approximated as a broken 
power law. We confirmed that this approximation 
is valid by using gtohssim to simulate a power law 
spectrum with photon-photon absorption above a 
threshold energy and showing that fitting a broken 
power law using gtlike correctly identified this en- 
ergy. Because gtlike can only fit a power law with 
a single break, to test the double-absorber model 
it was necessary to split the photon data from each 
object into low and high energy sets which were 
then each fitted with a broken power law. 

For the low energy set we calculated the mid- 
point between the two predicted breaks in the ob- 
server frame and took all photons below this en- 
ergy. This set was then fitted with a broken power 
law and tested against a simple power law null- 
hypothesis. The optimum break energy was then 
tested to see if it was consistent with the double- 
absorber model which predicts a break between 4 
and 7 GeV in this set. We then selected any pho- 
tons above the best fit break energy as a high en- 
ergy set. We fitted a broken power law to find any 
optimum second break in the spectrum and tested 
this against a simple power law null-hypothesis. 
Again, we also tested any break to see if it was 
consistent with the double-absorber model which 
predicts a second break between 19 and 30 GeV. 



els (jBurnham fc Andersonll2001l ). The AIC test is 
particularly suited to multi-model testing as it bal- 



^http: / /fermi. gsfc.nasa.gov/ssc/data/access/lat/BackgroundModels. html 
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2.4. Testing the Stability Of Spectra 

Assuming a break is present in tlie spectrum 
of an object, tiie stability of the break can give 
important information regarding its cause. In the 
double-absorber model the break energy is deter- 
mined by the wavelength of recombination lines 
and is therefore predicted to be stable. C onversely, 
in the model of Finke fc Dermer (I2OIOI) the break 
energy can vary on the timescale of several months 
as the emission region progresses along the jet. 
In general, variation in the spectral shape is ex- 
pected if it is determined by varying instantaneous 
physical parameters of the emission region (Klein- 
Nishina effects, B-field etc.). 

We split the dataset from each object into two 
equal time bins, epoch 1 and epoch 2, and found 
the spectral shape that best fitted each epoch. As 
a test of the stability we compared the break en- 
ergies of the broken power law fits for each object 
at the two epochs to see if they were consistent. 

3. Results 

For illustrative purposes, fits of the different 
spectral forms to each object are plotted along 
with aperture photometry data in Figure [TJ The 
aperture photometry data shows the differential 
photon flux in energy bins for each object. For 
each energy bin we took all photons reconstructed 
to originate within an on-region of 1° radius cen- 
tered on the source. We then calculated the 
flux using gtexposure to find the average expo- 
sure for each energy bin and subtracted an off- 
region to account for the background. The off- 
region was of the same radius as the on-region 
centered within 15° from the target and at least 
3° remov ed from any source s in the Fermi 1 Year 
Catalog (|Abdo et al.ll2010bl ). 



3.1. 



Testing For Deviation from a Simple 
Power Law 



The spectra of all of our sources show significant 
deviation from a simple power law. Tables [U [21 
and |3] show the parameters of the best fit to each 
source using a simple power law, broken power law 
and log-parabola fit respectively. 



3.2. Testing For Energy Breaks Caused By 
BLR Pair-Production 

Splitting the dataset for each object into a low 
energy and high energy set allowed us to test the 
predictions of the double-absorber model. After 
splitting the data in this way, the results for bro- 
ken power law fits to each source and the signifi- 
cance of the improvement they offer over a simple 
power law are shown in Tables U (low energy sets) 
and[5](high energy sets). Plots of the likelihood of 
a power law with a given break energy fitting the 
data are shown in Figures [2] (low energy sets) and[3] 
(high energy sets) . The regions where the double- 
absorber model predicts breaks to occur are shown 
as shaded boxes. 

When examining the low energy sets, in seven of 
the nine objects a broken power law is found to be 
a better description than a simple power law to > 
99% significance. However, in most of the objects, 
the break energy does not lie within the energy 
region predicted by the double-absorber model. 

Fits to the high energy set of most objects were 
not significantly improved by using a broken power 
law. Of the three objects that do reject the sim- 
ple power law to > 99% significance, the break 
energy is again not found in the region predicted 
by the double-absorber model. In PSIO evidence 
of absorption is found for seven objects in which 
we found a simple power law to be an adequate 
description, but the results of PSIO were limited 
by statistics. 

3.3. Testing the Stability Of Spectra 

Results from splitting the data set for each ob- 
ject into two epochs are shown in Table [6] Five of 
the nine sources were best fitted by different spec- 
tral shapes in each epoch (i.e. the best fit changed 
from a log-parabola to a broken power law or vice 
versa). For the broken power law fits, the break 
energies in six sources were found to differ over the 
two epochs by more than one standard deviation, 
and in three sources they were found to differ by 
more than two standard deviations. This is more 
deviation than expected from statistical variation 
alone, but with a sample size of only nine sources 
the evidence for variability is weak. 
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4. Discussion and Summary 

For each of the nine sources we took ah photons 
with energy below 12 GeV in the rest frame of 
the source and searched for a break. The double- 
absorber model predicts that a break should occur 
in this set between 4 and 7 GeV. Our results do 
not support this because although the fit to the 
spectra is in most cases significantly improved by 
introducing a break, the energy at which the break 
occurs does not fall in the predicted region. We 
then took all photons with energies above the low 
energy break and searched for a second break in 
the spectrum. The double-absorber model pre- 
dicts that such a break should occur between 19 
and 30 GeV. Again our results do not support 
this conclusion because fits to most objects were 
not significantly improved by introducing a break 
and in cases where they were the breaks did not 
fall within the predicted region. 

IStern fc Poutanen ( 2011 ) examined the pos- 
sibility of an intrinsic log-normal spectruno of 
7-rays that contains two 'breaks' due to absorp- 
tion with BLR recombination photons. As de- 
scribed in Section 12. 2[ we ran simulations of a 
log-parabola with a photon-photon absorption 
feature and found that fitting a broken power 
law would correctly identify the break ener- 
gies in a log-para b ola sp ectrum as proposed by 
Stern fc PoutanenI ( 2011 ). Our results do not 
support a double-absorbed log-parabolic spectrum 
since the predicted break energies are disfavored in 
the observational data. For BLR pair-production 
absorption to remain plausible the photoioniza- 
tion models used in PSIO would require significant 
modification. In several of the objects in our sam- 
ple, a break at ~ 1 - 3 GeV is seen. To create such 
a break, a region with a high ionization parameter 
line-like feature at '-^ 70 - 110 eV that dominates 
over H and He lines would be needed. 

When examining all of the data (not splitting 
the dataset by energy or time), the spectra of 
all the sources deviate significantly from a sim- 
ple power law. However, without using physi- 
cal motivation to prefer one model or the other, 
a broken power law is not always a better de- 
scription than a log-parabola. An AIC test indi- 
cates three objects, PKS 1502-^106, 3C 279, and 



RGB J0920-I-446, were described significantly bet- 
ter by a log-parabola than a broken power law, and 
PKS 0454-234 was described marginally better by 
a log-parabola. It might be expected that the log- 
parabola would be favored for weaker sources since 
with the normalization energy fixed as it is here 
it has one fewer degree of freedom than a broken 
power law. However; these four sources were not 
the faintest in our sample, in fact 3C 279 had the 
second greatest fiux. This suggests there is a gen- 
uine physical difference in these objects. 

For those objects that were better described by 
a broken power law, the break appears between 
~ 2 - 5 GeV in the source rest frame. The break 
energies determined in the low energy set and with 
no energy cut are consistent in all objects except 
3C 454.3, PKS 1502-^106, and RGB J0920-h446, 
which may suggest the presence of curvature in 
the spectra of these objects as well as any break. 
We found weak evidence that the break energies 
varied with time suggesting that the primary cause 
of the break is intrinsic to the emission region. 

Five of the nine sources are not best fitted by 
the same spectral shape in both time epochs (i.e. 
the fit changes from a log-parabola to a broken 
power law or vice versa). This could be a change 
in the spectral shape of the distribution of the par- 
ticles emitting the spectrum or it could be that 
when a spectrum with a given spectral shape and 
time-varying parameters is time-averaged it be- 
comes best fitted by a different spectral shape, a 
hypothesis that might be testable as more data 
become available. 

Both the wide range of break energies and 
the variability suggests that the cause of the 
break cannot be absorption with recombination 
line photons. lAbdo et al. ( 2010a ) note that the 
break energy is too low to be attributed to Klein- 
Nishina effects and conclude it is likely caused by 
a change in the energy distribution of electrons 
in the jet (assum i ng a le ptonic emission scenario). 
Finke fc Dermeil (|2010t ) claim this scenario pre- 



^They note that this is essentially the same function as a 
log-parabola in a different representation. 



diets a correlation between the spectral indices 
in the optical and 7-ray regimes, but this is not 
observed. They suggest instead two sources of 
seed photons for external Compton scattering can 
produce a sharp break, with Klein-Nishina effects 
providing an important contribution. In this sce- 
nario the break energy changes only gradually as 
the emission region progresses along the jet; how- 
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ever if a new component was ejected from the 
central black hole, an event thought to cause flar- 



ing in the radio and 7-ray regimes (jAcciari et al 



the break energy could be expected to alter 
rapidly. This view is consisten t with the observa- 
tions published bv lAbdo et all (I2OIII) . 

In summary, we find evidence that disfavors 
the double-absorber model for bright Fermi FS- 
RQs. The energy spectrum of 7-rays from bright 
FSRQ objects is sometimes better described by a 
log-parabola and sometimes by a power law with 
an energy break at a few GeV which varies with 
time. In order to investigate the causes of these 
features, future work is planned to extend the sam- 
ple size, to compare FSRQs to BL Lac objects and 
to investigate any correlations with flux. 
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Table 1 
Simple Power Law Fits 



\J DJCCt 


Index 

r 




3C 454.3 


2.485 ± 0.006 


653.5 


PKS 1502+106 


2.41 ± 0.02 


78.1 


3C 279 


2.41 ± 0.01 


70.8 


PKS 1510-08 


2.48 ± 0.01 


35.6 


3C 273 


2.76 ± 0.01 


32.9 


PKS 0454-234 


2.34 ± 0.02 


49.3 


PKS 2022-07 


2.52 ± 0.04 


7.8 


TXS 1520-1-319 


2.42 ± 0.01 


41.6 


RGB J0920-I-446 


2.35 ± 0.02 


43.1 



Note. — The best fitting parameters for simple 
power law fits to each object as described in Equa- 
tion AAICmin, is the difference in AIC between 
the simple power law fit to the object and the best 
fit out of simple power law, broken power law and 
log-parabola, a value of indicating a simple power 
law is the best fit to the data and a value > 2 indi- 
cating either a broken power law or a log-parabola 
is a significantly better fit to the data. 



Table 2 
Broken Power Law Fits 



Object 


Break Energy 


Indcxl 


Index 2 


AAIC„i„ 






Ti 


Ba 






(GeV) 








3C 454.3 




2.375 ± 0.007 


3.18 ± 0.03 





PKS 1502+106 




2.28 ± 0.02 


3.1 ± 0.1 


2.6 


3C 279 


9 1 +0.2 
^■-^-0.3 


2.32 ± 0.02 


2.70 ± 0.04 


6.6 


PKS 1510-08 


3-5+^5 


2.40 ± 0.02 


2.81 ± 0.06 





3C 273 




2.72 ± 0.02 


3.4 ± 0.1 





PKS 0454-234 


2 6+"-^ 

^■°-0.4 


2.12 ± 0.04 


2.61 ± 0.05 


1.8 


PKS 2022-07 


4 0+^ ** 


2.34 ± 0.07 


2.8 ± 0.1 





TXS 1520+319 


4 7+^-^ 
' -o.e 


2.36 ± 0.02 


2.76 ± 0.06 





RGB J0920+446 


16^5 


2.3 ± 0.1 


3.6 ± 0.4 


13.2 



Note. — The best fitting parameters for broken power law fits to each object as 
described in Equation[2] Sj, is in the object's rest frame. AAICmin is the difference 
in AIC between the broken power law fit to the object and the best fit out of simple 
power law, broken power law and log-parabola, a value of indicating a broken 
power law is the best fit to the data and a value > 2 indicating a log-paraboIa is a 
significantly better fit to the data. 



Table 3 
Log-Parabola Fits 



UDJCCt 


Index 

r 


Curvature 

p 




3C 454.3 


2.307 ± 0.009 


0.133 ± 0.006 


12.8 


PKS 1502+106 


2.02 ± 0.04 


0.16 ± 0.02 





3C 279 


2.36 ± 0.01 


0.071 ± 0.009 





PKS 1510-08 


2.35 ± 0.03 


0.07 ± 0.01 


2.4 


3C 273 


2.76 ± 0.02 


0.07 ± 0.01 


1.5 


PKS 0454-234 


2.11 ± 0.04 


0.12 ± 0.02 





PKS 2022-07 


2.27 ± 0.09 


0.11 ± 0.04 


0.4 


TXS 1520-1-319 


2.38 ± 0.01 


0.05 ± 0.01 


9.3 


RGB J0920-I-446 


2.26 ± 0.03 


0.10 ± 0.02 






Note. — The best fitting parameters for log-parabola fits to each 
object as described in Equation [3] /S.AlCmin is the difference in 
AIC between the log-parabola fit to the object and the best fit out 
of simple power law, broken power law and log-parabola, a value of 
indicating a log-parabola is the best fit to the data and a value 
> 2 indicating a broken power law is a significantly better fit to the 
data. 



Table 4 

Spectral Properties of Objects over the Whole Time Period at Low Energy 



Object 


Min Energy 


Peak log(iyF„) 




Index 1 


Index 2 


Sig over 


4.8 GeV 


4-7 GeV 






(GeV) 


(GeV) 


Pi 


r2 


simple PL 


exclusion 


exclusion 


30 454.3 


0.46 


0.244 


Q , + 0-2 


2.361+0.008 


2.86+0.03 


>99% 


>99% 


> 99% 


PKS 1502+106 


1.28 


0.934 


9 S + O-B 
^■ = -0.4 


2.16+0.06 


2.37+0.05 


90% 


86% 


93% 


30 279 


0.15 


0.101 


-^■"-0.4 


2.30+0.02 


2.60+0.04 


>99% 


>99% 


>99% 


PKS 1510-08 


0.41 


0.112 


4 0+°-'^ 

^■^-1.0 


2.41+0.02 


2.9+0.1 


>99% 


76% 


<1% 


30 273 


0.12 


0.005 


o -1+0.2 
''■-^-0.7 


2.70+0.02 


3.4+0.1 


>99% 


>99% 


>99% 


PKS 0454-234 


0.56 


0.523 


9 -1+0-3 
"'■-^-0.6 


2.13+0.04 


2.45+0.06 


>99% 


>99% 


>99% 


PKS 2022-07 


0.24 


0.693 




2.38+0.06 


3.1+0.4 


83% 


42% 


<1% 


TXS 1520-1-319 


0.29 


0.253 




2.36+0.02 


3.0+0.2 


>99% 


45% 


<1% 


RGB ,10920+446 


0.32 


0.139 


-1 O + 0.2 
-^■•'-0-2 


2.01+0.07 


2.42+0.04 


>99% 


>99% 


>99% 



Note. — Table of results over the whole available time period in the energy region from Min Energy to 12 GeV. Values for peak log{i>'Fiv) 
taken from lAbdo et al.l i2Q10cr except for PKS 2022-07 (see main text). Ehreaki is the value for the low energy spectral break which provided 
the best fit to the observed data. Fi and T2 are the photon spectral indices before and after the energy break respectively. Column 7 shows 
how significantly a broken power law improves the fit to the data compared with a simple power law. Columns 8 and 9 show the confidence 
level to which one can exclude energy breaks occurring at 4.8 GeV and in the 4—7 GeV break region as predicted in PSlO. All energies given 
in the object's rest frame. 



Table 5 

Spectral Properties of Objects over the Whole Time Period at High Energy 



Object 


Min Energy 




Index 3 


Index 4 


Sig over 


19.2 GeV 


19.2-30 GeV 




(GeV) 


(GeV) 


Ps 


r4 


simple PL 


exclusion 


exclusion 


3C 454.3 


3.3 


' ■•^-0.7 


2.69±0.06 


3.35±0.7 


>99% 


>99% 


>99% 


PKS 1502+106 


2.8 


12+t 


2.36±0.05 


3.1±0.1 


>99% 


94% 


>99% 


3C 279 


1.6 


35+1 


2.6±0.2 


>5 


98% 


94% 


95% 


PKS 1510-08 


4.0 


15if 


3.0±0.1 


2.5±0.2 


60% 


40% 


53% 


3C 273 


2.1 




3.3±0.2 


>5 


52% 


34% 


47% 


PKS 0454-234 


2.1 




2.50±0.07 


2.9±0.2 


87% 


93% 


99% 


PKS 2022-07 


6.0 




3.2±0.2 


1.7±0.5 


90% 


83% 


11% 


TXS 1520-1-319 


5.0 




2.7±0.1 


3.4±0.6 


47% 


69% 


63% 


RGB J0920-I-446 


1.3 


17tl 


2.39±0.04 


3.5±0.3 


>99% 


53% 


70% 



Note.— Table of results over the whole available time period in the energy region after the first break. Min Energy 
is the lowest photon energy for inclusion. Ebr&ak2 is the value for the high energy spectral break which provided the 
best fit to the observed data. Fs and T4, are the photon spectral indices before and after the energy break respectively. 
Column 6 shows how significantly a broken power law improves the fit to the data compared with a simple power law. 
Columns 7 and 8 show the confidence level to which one can exclude energy breaks occurring at 19.2 GeV and in the 
19.2 — 30 GeV break region as predicted in PSIO. All energies given in the object's rest frame. 



Table 6 

Spectral Properties of Objects over different time periods. 



Object 


Ebreak (GeV) 

No Time cut 


AAIC„i„ 


Eweak (GeV) 
Epoch 1 


AAIC„i„ 


Eweak (GeV) 
Epoch 2 


AAIC,„j„ 


3C 454.3 


4 5+"^^ 










4 5+°-^ 
^■°-o.i 


0.32 


PKS 1502+106 


11 2 + 1-1 

-1 -1-^-0. 8 


2.6 


10 7+1-* 

' —0.8 


0.5 


9 O + 0.4 
^■•^-0.3 





3C 729 


9 1+0.2 
^■1-0.3 


6.6 


1 7+0.5 
^- ' -0.1 


6.1 


2 4+°-<' 


2.9 


PKS 1510-08 


n c + 0.6 





o -1+0.8 
•^•-1-0.3 





3 9+°-^ 





3C 273 


9 9 + 0-2 
^■^-0-5 





J 3 + 0.4 
-^•■'-0.2 


8.9 


4 1+0-5 

^•-^-0.5 





PKS 0454-234 


2 6+°^^ 
^-°-0.4 


1.8 


o q + 0.3 
■^•■^-0.5 


1.5 







PKS 2022-07 







3 R + l 1 










TXS 1520-1-319 









0.7 


4 4+2.0 





RGB J0920-I-446 


16+4 





1 0+°* 
^•"-0.1 





r O+0.6 
0.0_j 2 






Note. — Table of results found by splitting the data for each object into two equal time epochs. The most likely break energy 
for each object, in the object's rest frame, with no time cut and in each epoch individually is shown along with 68% confidence 
intervals and the difference from the AIC minimum, a value of indicating a broken power law is the best fit to the data and a 
value > 2 indicating a log-parabola is a significantly better fit to the data. 



3C 454.3 PKS 1502+106 3C 279 




log|,|(Energy(GeV)) ]ogn,IEiiergy(GeV)) logmCEnei-gj-fGeV)) 



PKS 2022-07 TXS 1520+319 RGB J0920+446 




logL,i(Energy(GeV)) ]ogiii(Energy(GeV]) iogni(Energy(GeV)) 



Fig. 1. — Aperture photometry data for each source along with the following spectral shapes from a maxi- 
mum likelihood fit to the full dataset: simple power law (dotted), broken power law (dashed), log-parabola 
(solid), and double-a,bsorbe r (dot-dashed). Error bars are 68% confidence intervals for each model (see 
Aggarwal fc Caldwelll (|2012l )). Residuals ((Observed Flux - Model Flux) /Observed Flux) for the best fitting 



spectral shape are shown underneath each plot. Energies are in the observer frame. 
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Fig. 2. — Plots of -log(Likelihood) against power law break energy, in the object's rest frame, for the low 
energy datasct of each object. The shaded box shows the region in which the double-absorber model predicts 
an energy break to occur. 
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5 10 15 20 25 30 
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PKS 1510-08 
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40 
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Energy (GeV) 





10 15 20 25 30 
Energy (GeV) 



T3.91993e5 



PKS 0454-234 
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Energy (GeV) 



RGB J0920+446 




10 20 30 40 50 60 
Energy (GeV) 



Fig. 3. — Plots of -log(Likelihood) against power law break energy, in the object's rest frame, for the high 
energy datasct of each object. The shaded box shows the region in which the double-absorber model predicts 
an energy break to occur. 
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